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(N ■ Abstract 



We consider a queuing model with the workload evolving between consecutive i.i.d. exponen- 
tial timers {e^ }j=i,2,... according to a spectrally positive Levy process Y(t) which is reflected 

(i) 

at 0. When the exponential clock e q ends, the additional state-dependent service require- 
ment modifies the workload so that the latter is equal to Fi(Y(eq)) at epoch + ■ ■ ■ + dp 
■ for some random nonnegative i.i.d. functionals Fi. In particular, we focus on the case when 

p j Fi(y) = (Bi — y) + , where {-Bi}j=i,2,... ar e i.i.d. nonnegative random variables. We analyse the 

^ 1 steady-state workload distribution for this model. 

> 



1 Introduction 

In this paper we focus on a particular queuing system with additional state-dependent services. 
There has been considerable previous work on queues with state-dependent service and arrival 
processes; see for example the survey by Dshalalow [12] for several references. The model under 
consideration involves a reflected Levy process connected to the evolution of the workload. Special 
cases of Levy processes are the compound Poisson process, the Brownian motion, linear drift 
£Sj ■ processes, and independent sums of the above. The literature on queueing systems driven by Levy 

! processes is rather limited; see e.g. Bekker et al. [3] for references. 

Specifically, in this paper we consider a storage/workload model in which the workload evolves 
according to a reflected at zero spectrally positive Levy process Y(t). That is, let X(t) be a 
spectrally positive Levy process (a Levy process with only positive jumps) modelling the input 
minus the output of the process and define — inf s ^o X(s) = and X(0) = x ^ 0. Then we 
have that Y(t) = X(t) - wi ss gX(s) (where Y(0) = x for some initial workload x ^ 0). In 
addition, at exponential times with intensity q, given by {eq ^}i=i : 2,..., the workload is "reset" to a 
certain level, depending on the workload level before the exponential clock ends. Specifically, at 
epoch t = e q ^ + • • • + eq the workload V(t) equals Fi(V(t— )) for some random nonnegative i.i.d. 
functionals i^. 

The main goal of our paper is to derive the stationary distribution of the workload V{t) for 
the above-described queuing model. We first identify the stationary distribution of the workload 
at embedded exponential epochs and then extend this result to an arbitrary time by using renewal 
arguments. We also identify the tail behaviour of the steady-state workload. 
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This model unifies and extends several related models in various directions. First of all, if X 
is a compound Poisson process and if F{ is the identity function, then our model reduces to the 
workload process of the M/G/l queue. Kella et al. [18] consider a model with workload removal, 
which fits into our model by taking Fi(x) = and by letting the spectrally positive Levy process 
X be a Brownian motion superposed with an independent compound Poisson component. The 
added generality allows one to analyse more elaborate mechanisms of workload control, where the 
exponential times can be seen as review times, during which the workload can be changed to a 
different level as desired. Allowing for general functions Fi opens up the possibility of optimising 
such controls, although we do not consider this problem here. Instead of considering the classical 
compound Poisson input process and a linear output process for the queuing model, one can consider 
the case in which the Levy process is nondecreasing, i.e. a subordinator; see for example Bekker et 
al. [3] and Boxma et al. [7]. 

Apart from the wish to unify and extend clearing models, we were also challenged by introducing 
a continuous-time analogue of the alternating service model considered in Vlasiou et al. [28, 29, 30, 

31, 33] that gave rise to the Lindley-type recursion W = (B - A — W) + . In particular, we focus on 
the case when F{(y) = (Bi — y) + , where {-Bj}i=i,2,... are i.i.d. nonnegative random variables. Hence, 
at exponential epochs the controlling mechanism leaves only a portion of the workload depending 
on the size of the workload just prior to the exponential timer. In particular, if Fi(y) = (Bi — y) + , 
then this mechanism keeps the workload below a generic random size B, decreasing it when it is 
relatively large at the exponential epoch and increasing it when it is much smaller than B. This 
can be viewed as a continuous-time analogue of the above mentioned Lindley-type recursion. In 
the context of workload control mentioned above, this example can be interpreted as a mechanism 
that reflects existing storage with respect to an upper bound B. 

Our results focus on qualitative and quantitative properties of the steady-state workload distri- 
bution. We first establish Harris recurrence for the Markov chain embedded at workload adjustment 
points, yielding the convergence of the workload processes to an invariant distribution. We derive 
an equation for the invariant distribution of the embedded chain, as well as the invariant distribu- 
tion of the original process. We use this equation to obtain expressions for the invariant distribution 
for an example that generalises [3, 7, 18]. We also investigate the tail behaviour of the steady-state 
distribution under both light-tailed and heavy-tailed assumptions. All these groups of results have 
the common theme that we rely on recently obtained results in the fluctuation theory of spectrally 
positive Levy processes. 

The paper is organised as follows. In Section 2 we introduce a few basic facts concerning 
spectrally positive Levy processes. In Section 3 we consider the embedded workload process and 
derive a recursive equation for its stationary distribution. In Section 4 we determine the steady- 
state workload distribution. Later on, in Section 5 we present some special cases. Finally, in Section 
6 we focus on the tail behaviour of the steady-state workload. 

2 Preliminaries 

Throughout this paper we exclude the case of X with monotone paths. Let the dual process of 
X(t) be given by X(t) = —X(t). The process {X(s),s ^ t} is a spectrally negative Levy process 
and has the same law as the time-reversed process {X((t — s)—) — X(t), s ^ t}. Following standard 
conventions, let X(t) = inf s ^A(s), X(t) = sup s < t X(s) and similarly X(t) = ini s <^t X (s), and 

X(t) = sup s ^ t X(s). One can readily see that the processes Y(t) = X(t) — X(t) (for Y(0) = 0) 
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and X(t) (where X(0) = 0) have the same distribution; see e.g. Kyprianou [20, Lemma 3.5, p. 74]. 
Moreover, 

-X(t) = X(t), X(t) = -X(t). 

Since the jumps of X are all non-positive, the moment generating function E[e sx ^] exists for 
all 9 and is given by E[e ex V] = for some function ip(9) that is well defined at least on the 

positive half-axis where it is strictly convex with the property that lim^oo i^(0) = +oo. Moreover, 
ip is strictly increasing on [$(0),oo), where <£(0) is the largest root of ip{0) = 0. We shall denote 
the right-inverse function of ip by $ : [0, oo) — > [^(0), oo). 

Denote by a the Gaussian coefficient and by v the Levy measure of X (note that a is also a 
Gaussian coefficient of X and that n(A) = v{— A) is a jump measure of X). Throughout this paper 
we assume that the following (regularity) condition is satisfied: 

,o 

a > or / xu(dx) = oo or u(dx) << dx, (2-1) 

where « dx means absolutely continuity with respect to the Lebesgue measure. Moreover, we 
assume that 

P,(r " < oo) = 1, (2.2) 

where 

Tq = M{t > : X(t) s£ 0}. 

Finally, P x denotes the probability measure P under the condition that X(0) = x, and E x indicates 
the expectation with respect to P x . 



2.1 Scale functions 

For q 0, there exists a function : [0, oo) 

continuous and increasing with Laplace transform 



[0, oo), called the q-scale function, that is 



poo 

/ e- e vW®(y)dy = M0)-q)-\ (2-3) 
J o 

The domain of is extended to the entire real axis by setting W^ q \y) = for y < 0. We 

mention here some properties of the function that have been obtained in the literature which 
we will need later on. 

On (0, oo) the function y i— > W^ q \y) is right- and left-differentiable and, as shown in [23], under 
the condition (2.1), it holds that y i— > W^ q \y) is continuously differentiable for y > 0. 
Closely related to is the function Z^ 9 ' given by 

Z<s\y) = l + q I 'w iq \z)dz. 
Jo 

The name "g-scale function" for W^ 9 ' and is justified as these functions are harmonic for the 
process X killed upon entering (— oo, 0). Here we give a few examples of scale functions. For a large 
number of examples of scale functions see e.g. Chaumont et al. [10], Hubalek and Kyprianou [17], 
Kyprianou and Rivero [22]. 
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Example 1. If X(t) = o~B(t) — fit is a Brownian motion with drift fx (a standard model for small 
service requirements) then 



where S = a 2 \J fx 2 + 2qo~ 2 and ui = fx/a 2 . 
Example 2. Suppose 

N(t) 

where p is the speed of the server and {<7j} are i.i.d. service times that are coming according to the 
Poisson process N(t) with intensity A. We assume that all o~i are exponentially distributed with 
mean 1/ fx. Then if) (6) = p9 — X8/(fi + 6) and the scale function of the dual W^ q > is given by 

W®{x) = p~ x (A+e q+ ^ x - A_e q ~^ , 
where A± = q +^_^^ with q + (q) = &(q) and q~{q) is the smallest root of ip(6) = q: 



± q + \-fip±y / (q + \- ftp) 2 + Apqfx 
q {q) = 2p • 

2.2 Fluctuation identities 

The functions W^ q ' and play a key role in the fluctuation theory of reflected processes as shown 
by the following identity (see Bertoin [4, Theorem VII. 4 on p. 191 and (3) on p. 192] or Kyprianou 
and Palmowski [21, Theorem 5]). 

Lemma 2.1. For a > 0, 

E ( e - a X(e q )\ = g(«-$(g)) 



*(g) & (a)-qY 



which is equivalent to 



P{X{e q ) € dx) = ^J^ wiq) ( dx ) ~ qW {q \x)dx, x > 0. 

Moreover, —X_(e q ) follows an exponential distribution with parameter $(<?). 

The scale function gives also the density r^ q \x,y) = R^ q \x,dy)/dy of the g-potential measure 



POO 

R®(x,dy) := / e-*P x {X(t) £ dy,r > t)dt 
Jo 



of the process X killed on exiting [0, oo) when initiated from x. See also Pistorius [26]. 
Lemma 2.2. Under (2.1), we have that 



-fe)(x,i/)= / e~*M z W {q)l {y-x + z)-<S>{q)W {q \y — x + z) 

J[(x-y)+,x] 



dz. 
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Proof. We start by noting that for all x, y > and q > 0, 



R^(x,dy) = ^P x (X(e g ) G dy,X(e g ) > 0). 

From the Wiener-Hopf factorisation of the Levy process we have that X{e q ) — X_{e q ) is independent 
of X(e q ). This leads to 

Rb)(x,dy) = -P{X(e q ) G dy - x,X(e q ) ^ -x) 

= -P({X(e q ) - X(e q )) + X(e q ) G dy - x, -X(e q ) ^ x) 



- / P(-K(e q ) G dz)P(X(e q ) - X(e q ) edy-x + z). 

Q J\(x-v)+,x] 



Q J[(x-y) + ,x] 

In the above expression, we integrate over the position of — X_(e q ) which is not less than under 
P x (this leads to the condition that — X_{e q ) $5 x under P = Pq) and it is less than X(e q ) = y under 
P x (hence, — X_(e q ) > x — y under P). Note that we always have that — X(e q ) ^ under P, and 
thus the above integral is equal to the integral over [0, x] when y > x. 

Recall, that by duality X(e q ) — X_(e q ) is equal in distribution to X(e q ) which has been identified 
in Lemma 2.1. In addition, the law of — X_(e q ) is exponentially distributed with parameter Q(q). 
We may, therefore, rewrite the expression for PS q >(x,dy) as follows: 

R^(x,dy) = I e-^ q)z \w {q) (dy-x + z) - <£>(q)W {q) (y - x + z)dy\ dz. (2.4) 

J[(x-y)+,x] 

Under Condition (2.1), W^ q ' is differentiable and hence the last equality completes the proof. □ 

Remark. Lemma (2.2) and similar results can be proven without the assumption made in (2.1), 
but at the cost of more complex expressions. We would have to use (2.4) instead of the much nicer 
form for A q \x,y)dy. 

3 Equilibrium distribution of the embedded process 

We consider the workload process at the embedded epochs e q ^ + • • • + e q n \ just after the additional 
service arrives. Note that this process is a Markov chain {Z n ,n G N} with transition kernel: 

k(x,dy) = P x {F{Y{e q )) G dy) = j P x (f(Y(e q )) G dy)dP F (f), (3.1) 

where P F is the law of F. 

Lemma 3.1. We have that P x (Y(e q ) G dy) = h(x,y)dy + e~^ q >W^ (0)8 (dy) , where 



h(x,y)=qr^(x,y)+e-^ q)x 



.*(?) 

and where the first increment qA q \x,y) is given in Lemma 2.2. 



W®'(y)-qW®(y) 



(3.2) 
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Proof. Define kq = inf{t ^ : Y(t) = 0}, and observe that 

P x {Y(e q ) G dy) = P x (Y(e g ) G dy, k > e q ) + P x (Y(e q ) G dy, k < e q ) 

= P x {X{e q ) G dy,T - > e ? ) + P(F(e g ) G dy)P x {r~ < e q ) 
= qrM(x, y)dy + P(X(e q ) G dy)P(-X(e ? ) > x) 



qr^(x,y)dy + e-^ x 
+ e-*W x wM(0)6 (dy) 



!{y>o}dy 



where in the second equality we use the lack of memory of the exponential distribution and the 
fact that X is spectrally positive; hence it crosses in continuous way. The last equality follows 
from Lemma 2.1. □ 

Lemma 3.2. Assume that there exists an a.s. finite r.v. Fq such that 

F(y) < F for any y. (3.3) 

Then the stationary distribution ir(y) of Z n exists and satisfies the following balance equation: 



g{y)dn(y) + g(0)7r(0) = / / g{y)k(x, dy)dvr(x) (3.4) 

J[0,oo) J[0,oo) 

for any bounded function g. 

Proof. We show that {Z n ,n G N} is Harris ergodic, see for example Asmussen [2, Theorems 
VII.3.3 and VII.3.5, p. 200-201]. Fix N > such that P(F < N) > and define t n = inf{n ^ 

1 : Z n ^ N}. It is easy to construct i.i.d. random variables {Fo tn ,n G N} such that p^i = Fq and 
Z n ^ Po,n> G N. Observe that 

P(t n >k\Z = x)= P(Z l >N,...,Z k >N\Z = x)^ P(F 0A >N,..., F , fc > N) = P(F > N) k , 
which implies that 

sup E[tn I Z = x] < oo. (3-5) 

This implies Harris ergodicity, once we find a constant p > and a probability measure Q(-) such 
that 

P(Z 1 £ B \ Z = x) ^ pQ{B), x€[0,N]. (3.6) 

Let x G [0, iV]. We construct p and Q(-) as follows. Recall that kq = inf{i ^ : = 0}, observe 
that 

P(Zi G P | Z = x) = P c (P(F(e ? )) G P) 

>P x (F(Y(e q ))eB, Ko <e q ) 
= P(F(Y(e q ))eB)P x (K <e q ) 
> P(F{Y(e q )) G B)P n (k < e q ) 
■= Q(B)p. 

Since the paths of Y(-) are non-monotone, we have that p > 0, implying (3.6). 

□ 
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Remark. The assumption (3.3) is satisfied when F(y) = (B — y) + with Fq = B. Moreover, if the 
functional is given by F{y) = (B — y) + and the r.v. B has a density then we have that 



7T(0) 



P f'OO f'OG f 

/ / / \qr ( - q \x,y)dy + e 

J\0,oo)J0 Jt I 



'[0,oo) 

and for any bounded function g, 



W®'(y)-qW®(y) 



dy}dF B (t)dir(x), 

(3.7) 



g(y)dir(y) 

roo ft 

[0,00) y Jo 



dyg(t-y)\qrM(x,y)+e-^ x 



$(<?) 
+ W^(0) 



W {q) '{y) -qW {q \y) 



dF B (i)dvr(x) 



[0,00) 



*W*d7r(x) / 5 (t)dF B (t), 



where .Fg is the distribution function of the generic r.v. B. 
In Section 5 we analyse more specific examples. 



4 Steady-state workload distribution 

Theorem 4.1. Suppose that the stationary distribution ir exists. Then the stationary distribution 
V(°°) a ^ so exists. Moreover, for a bounded function g, 



Eg(V(oo)) 



[0,00) 



7r(dx) / g{y)dy{ qA q) (x,y) +e 



-$(g)a; 



W (<?) '(y) -^ (9) (y) 



+ g(P)wM(0) 



n(dx), 



where the distribution ir satisfies (3.4) with k defined via (3.1) and Lemma 3.1 and where A q \x,y) 
is given in Lemma 2.2. 

Proof. The first part of the theorem follows from Asmussen [2, Theorem VII. 6. 4, p. 216]. Using the 
Palm inversion formula (see Asmussen [2, Theorem VII. 6. 4, p. 216]) one can identify the steady- 
state workload distribution V(oo) by the following identity: 



Eg{V{oo)) = q / n{dx)E x \ g(V(s))ds. 
J [0,00) 



(4.1) 
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The RHS can be further developed as follows: 

q [ n(dx)E x r g(V(s))ds = q [ n(dx)E H E x [g(Y(s))]ds 

J[0,oo) jO J[0,oo) io 

= q / ir(dx) / e- qt E x [g(Y(tWt 

J\0,oo) JO 



[0,oo) 

n(dx)E x [g(Y(e q ))} 



[0,00) 



7r(dx) / g(y)dy I qr 



'[0,oo) JO 

+ 5(0)^(0) / e-^ q)x 7r(dx). 



*(?) 



00 







□ 



5 Computational examples 

We now turn to analysing a few specific examples. We find that there are several solution strategies. 
One can either solve the equations given in Sections 3 and 4 directly, or one can also take a less 
direct route, using Laplace transforms. We shall consider examples of both strategies. 

To this end, we start with the following simple, but very useful observation. Using PASTA, V has 
the same distribution as U which is the equilibrium distribution of the Markov chain {U n ,n € N} 
of the workload process embedded at times (e^ + • • • + e^)— (i.e. right before the "correction"). 
To obtain our main result, we first state and prove the following useful lemma. 

Lemma 5.1. The following equality holds in distribution: 

U = max{F(U)+X(e q ),X(e q )}. (5.1) 
Proof. If Uq is x, we see that U\ = Y(e q ), with V(0) = F(x). If V(0) = F(x), we further have 

Y(t) = max{F(x) +X{t), X(t)}, 

which follows, for example, by mimicking the proof of Kyprianou [20, Lemma 3.5, p. 74], starting 
from the expression for Y(t) given in Kyprianou [20, p. 19]. Combining these two observations 
leads to the statement of the lemma. □ 

Example 3. The most trivial example is when F(y) = B ^ 0. In this simple case, there is no need 
to use the formula derived for the generator k(x,y). Using the above lemma, we see that 

V(oo) = max{£ + X(e q ),X(e q )} 

= X(e q ) + max{5 + X{e q ) - X(e q ), 0} 
v 



X{e q ) + max{S - e$( g ) , 0}. 



S 



In the last equation, which follows from the Wiener-Hopf factorisation, e$f q \ is a random variable 
which is exponentially distributed with rate &(q), which is independent of everything else. Observe 
that 



E j e - S max{B- e$(9) ,o}] = P fa {q) > B ) + £[e" s ( B - e *M > ; e* ( , } < 5] 

= P(e* ( ,) > B) + E[e-< B - e *(^} - E[ e -< B ~ e ^;e m > B] 

_ $(g)^[e~ 5B ] - s£[e-' I> fa) g ] 
*(?) - « 

Combining this with Lemma 2.1, we obtain 

Me -.V(oo)i = Q(s-^(q)) <S>( q )E[e-° B ]-sE[e-^ B } 

1 1 (*)-?) 



(5.2) 



This is an extension of various results in the literature focusing on clearing models (i.e. systems 
with workload removal) where B = 0. See for example [8, 18] and references therein. 

Example 4. We now consider an example where it seems more natural to solve the equations 
developed in Sections 3 and 4 directly. Consider the case when F{y) = (B — y) + with B being 

exponentially distributed with intensity (3. Moreover, X(t) = X^i*^ °i — pt is a, compound Poisson 
process with exponentially distributed service times Oi with intensity \i (see also the setup of 
Example 2). Note that 

/>oo 

$(0) = P 9-X (1 - e - 0z )/ie-^dz = P 6-\ 
Jo 

and recall that &(q) = q + (q)- Thus, 



Ee-^ = g [ .(dx) r d ye-^\ X ,y) + * ( «» + V + - A-M*®) 

J[o,oo) Jo ^(?)(^(a) -<?) p 



'[0,oo) 

H{a,n) + r—^r- — + -{A + - A-)ir{9{q)), 



where 



TT (a \ g +/ q + {q)-q (q) ~, Q , 2q + (q) 

H(6,u) = -A <u(q + {q))— ITTTT — it(0) 



p I V " ™"(e-q+(q))(e -q-(q)) W 2 -<?+(g) 2 



+^ + g + ((?)-g-(g)) 



2 -g-(g) 2 



and u(0) = Jjo e~ u(dx). To complete the computations we have to find the LST ff of the 
stationary distribution 7r. By the memoryless property of the exponential distribution of B we 
have that 7r(dx) = /3e~^ x dx, x ^ 0. Hence, 

7f(0) =7T(0) + 
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We will find now tt(0) using (3.7). 



tt(0) = tt(0)/3 



+ /3 2 



.^(<?)'(y) _^(?)(y) 



dy+ 



•^dt / <! gr( ? )(x,i/)+e-^ 



Let now e j a(dx) = /?e @ x dx and 



G(g) 



*(?) 

q~(q) - q + (q) 



W®'(y)-qW®(y) 



dy. 



p(P-q (<?)) 9 (q)q + {q) 



Then, 



vr(O) 



G(<l)v&®+H(P,ep)-H(J3,ep) 



1 - G(g) 

More general cases can be handled at the cost of more cumbersome computations. For example, we 
can add a compound Poisson process with phase-type jumps to the Levy process, and we can allow 
B to have a phase-type distribution. See [9] for similar computations in a discrete-time setting. 

If one cannot expect to obtain distributions in closed form, one can still aim to obtain Laplace 
transforms. Using similar arguments as in Example 3, we obtain the key equation (abbreviating 
V = V(oo)) 



E[e~ sV ] 



q{s - $ (g)) $(q)E[e- 8F W] - sE[e-*M F W] 



(5.3) 



$ (?) O («) - q) *(?) - s 

This equation is, of course, too complicated to solve for an arbitrary F, but nevertheless seems 
useful. 

Example 5. Suppose that F(x) = 5x, 5 G (0,1). This case is a generalisation of a model for 
the throughput behaviour of a data connection under the Transmission Control Protocol (TCP) 
where typically the Levy process is a simple deterministic drift; see for example [1, 16, 24, 25] and 
references therein. 

Equation (5.3) reduces to 



E\c 



-sVi 



E[e- sSV ] + 



qs 



Hq)W(s)-q) 



E[e~^ sv ]. 



(5.4) 



This is an equation of the form v(s) = g(s)v(5s) + h(s), which, since v(0) = 1, has as (formal) 
solution 



k-l 



v{s) = n 9(S j s) + ]T h(S k s) J] g(Ps) 

3=0 k=0 j=0 

Specialising to our situation we obtain 



v(s) 



n 

3=0 
oo 

n 

j=0 



— )+ vmq))E 



qs5 k 



k-l 



^to fjo^-^^y 
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Since ip(0) = 0, it easily follows that the infinite products converge, and the final expression for 
v(s) yields an equation from which the only remaining unknown constant v(5&(q)) can be solved 
explicitly 



6 Tail behaviour 

In this section we consider the tail behaviour of V = V(oo), under a variety of assumptions on 
the tail behaviour of the the Levy measure v. The treatment in this section can be seen as the 
continuous-time analogue of the results in Vlasiou and Palmowski [32]. There, a similar result 
is shown where e q is geometrically distributed, X(-) is replaced by a (general) random walk, and 
F(y) = {B — y) + , with B identical to the increments of the random walk. In [32] we modify 
ideas from Goldie [15] in the light-tailed case and develop stochastic lower and upper bounds in 
the heavy-tailed case. Here we take a different approach which is based on the well-developed 
fluctuation theory of spectrally one-sided Levy processes. Before we present our main results, we 
first state some lemmas. 

Lemma 6.1. The following (in) equalities hold: 
P(U >x ) = P(X(e q ) + F(U) > x) 

roo 

+ / (P(X(e g ) >x)- P(X(e q ) >x + y))P(-X(e q ) - F(U) e dy), (6.1) 
J o 

P(U > x) sC P{X{e q ) + F(U) >x)+ P(X(e q ) > x)P{-X{e q ) > F(U)), (6.2) 

P(U >x)> P{X{e q ) > x)P(-X(e q ) > F(U)). (6.3) 

Proof. All identities follow from Lemma 5.1, the decomposition X[e q ) = X(e q ) — (X(e q ) — X(e q )), 
and the Wiener-Hopf factorisation factorisation which implies that X{e q ) and X{e q ) — X(e q ) are 
independent. We also use the simple fact that X(e q ) — X{e q ) has the same law as —X_{e q ) (see 
Kyprianou [20, p. 158]). □ 

In addition, we need two standard results from the literature on Levy processes. 

Lemma 6.2. (Kyprianou [20, p. 165]) The random variable X(e q ) has the same law as H(e K ^ q ^), 
where {(L~ 1 (t), H(t)),t ^ 0} is a ladder height process of X with the Laplace exponent k(q,C) 
defined by 

E ^L- 1 (t)+CH(t) = e «(e,C)*_ 

Prom the above, one can easily derive the following version of the Pollaczek-Khinchine formula. 
Lemma 6.3. (Bertoin [4, p- 172]) The following identity holds: 

P(X(e q ) >x) = K(q,0)U iq) (x,oo), 

where 

/*oo roo 

U^{dx)= / e~ qs P{H{t) € dx,L~ 1 (t) G ds)di 
Jo Jo 

is the renewal function of the ladder height process {(L -1 (t) , H (t)) , t < L(e q )} and L(-) is a local 
time of X. 
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We now turn to the tail behaviour of U. Let 11(A) = v(—A) be the Levy measure of the 
spectrally positive Levy process X (with support on R+). First we investigate the case where the 
Levy measure is a member of the so-called convolution equivalent class S^ a \ To define this class, 
take a ^ 0. We shall say that measure LT is convolution equivalent (LI 6 S^ a ') if for fixed y we have 
that 

n(u — y) 

lim — =— - — = e ay , if LT is nonlattice, 

u->co H(u) 

n(n — 1) 

lim — - — = e a , if LT is lattice with span 1, 

LI(n) 



n— »oo 



and 



lim = 2 / e^n(dy), 

t^oo U(U) Jo 



where * denotes convolution and n(-u) = II((u, oo)). When a = 0, then we are in the subclass of 
sub exponential measures and there is no need to distinguish between the lattice and non-lattice 
cases (see [6]). We start from the following auxiliary result, which is the continuous-time analogue 
of Lemma 2 in [32]. 

Lemma 6.4. Assume that LT £ and ip(a) < q for ip(a) = \ogEe aX ^\ Then 

P{X(e q )>x)~ (g _' (Q))2 n(x), (6.4) 



g(g)_+ «i 
(g-V(a)) 2 ' <%) 



P(X( e? ) > x) ~ 7^-^5(1), (6.5) 



where /(x) ~ g(x) means that ]im. x ^ 00 f (x) / g(x) = 1. 
Remark. Note that for a = 

P(X(e 9 ) > x) ~ P(X(e 9 ) > x) ~ ^IT(x). 

Proof. It is well known that P(X(i) > x) ~ i£e a ^wrjjf(i)(x) for i fixed as x — > oo, where ^-x{t) 
is a Levy measure of X{t) (see Embrechts et al. [14]). Since X(t) is infinitely divisible we have 
Hx(t)(0 = *n(-) and hence P(X(i) > x) ~ t( J Be aX ( 1 )) t fl(x). Since X(e q ) < X(e g ) by (6.5) and the 
dominated convergence theorem we obtain (6.4). We will use similar arguments as in the proof of 
Lemma 3.5 of Kliippelberg et al. [19]. For LI# S note that 

P{H{t) > u) ~ tOEe^W/Stftu), 

where II# is the Levy measure of the process {H(t),t < e K ( q \} (see Embrechts et al. [14]). Using 
uniform in u Kesten bounds [19]: 

P{H{t) >u)^ P(H([t] + 1) > u) < K(e)(Ee aH W + e^n^u) 

for any e > and some constant K(e), and the dominated convergence theorem, we derive by 
Lemma 6.2, 

P(X(e q ) > u) 40) 

n H («) (4,o) -log ^e-^ 1 )) 2 ' y ' 
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The Wiener-Hopf factorisation yields that 

where (L _1 (t), H{t)) is a downward ladder height process with Laplace exponent k(g, £). Since X 
is spectrally positive, we can choose the process H (t) = —t and hence 

k(q,0) 



poo 

Ee aH( - e *^ = k(q,0) / e _/i(9 ' 0) *e- at dt 
Jo 

where k(q, a) = Q(q) + a. Thus 

q k(q,a) «(g,0) 



9 - «(?, 0) «(?, 0) - log Ee aH W ' 

Using the well known fact that q = K(q,0)k(q,0) (see Kyprianou [20, p. 166]) we identify the 
right-hand side of (6.6) as 



lim 



P(X(e q )>u) _ q (*(?)+ a) 



U H (u) (q-Ha)) 2 HQ) 

Now using similar arguments like in Vigon [27] (see also Kyprianou [20, Th. 7.7 on p. 191] and 
Kyprianou [20, Th. 7.8 on p. 195]) we derive 



poo 

Il H (u) = / U(u + y)V(dy), 
Jo 



where V(y) is the renewal function of the downward ladder height process {(L 1 (t),H(t)),t < 
k(q,0)} = {(L" 1 ^), H(t)),L~ l (t) < e q ). Thus 

«->oo U{U) Jo 

Ee - g L-Ht)-aH(t) dt 

1 1 



e -k(q,a)t dt 



k(q,a) $(q) + a' 

Hence, by [13] also Uh G if and only if LT £ . This completes the proof. □ 

It is known [14] that if for independent random variables \i (* = 1>2) we have P(xi > u ) ~ 
CiG(u) as u — > oo and G & S^ a \ then P(xi + X2 > u) ~ {c\Ee aX2 +C2Ee axi )G(u). This observation 
and (6.1) in Lemma 6.1 and Lemma 6.4 yield the following main result. 

Theorem 6.1. Assume that U £ and ip(a) < q. Moreover, let F{y) ^ Fq(^ 0) for any y, and 
assume that there exists a constant c such that P(F(y) > x) ~ P (i*b > x) ~ cll(a;) as x — > oo 
/or eac/i y (If c = then P{F{y) > x) = o(YL(x)) ). Then 

P{U >x)~ IcEe *^ + t ?-—^Ee aF W 

I (g-w*)r 

e -a(-X(e 9 )-F(t/);A . _x(e 9 ) - F(17) > ol j fl(x) 



g $(g)+a 
+ (g-^(a))2 



as x — > oo. 
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The conditions in this theorem are satisfied by both examples F(y) = (in which case we take 
Fq = 0, c = 0) and F(y) = (B — y) + (in which case Fq = B). If II is subexponential (IT € S^), 
then 

1' 



P(U >s)~(c + -l n(x). 

We will consider now the Cramer case (light-tailed case). Assume that there exists $((?) such 
that 

^{q))=q (6.7) 

and that 



m(q) :- 



OP 



< oo. (6. 

/?=-<%) 



Note that if II G S(a) and ip(ce) < q, then condition (6.7) is not satisfied. Moreover, we assume 
that 

Ee mnu) < (6 .9) 

Theorem 6.2. Assume that (6.7)-(6.9) hold and that the support ofU is non-lattice. Then 

P{U > x) ~ Ce-*ti* x 

as x — > oo, where 

C = P(-X(e q ) > F(U))K(q,0) (^(g)m(g))" 1 . 
Proof. We introduce the new probability measure 

dP e 



dP 



e ex{t)-^{e)t 



where Tt is a natural filtration of X. On P e , the process X is again a spectrally positive Levy 
process with the Levy measure Il^dx) = e dx U(dx), which is also nonlattice. Let be the 
renewal function appearing in Lemma 6.3 with P replaced by P 9 . Recall that L~ l {t) is a stopping 
time. Hence, from the optional stopping theorem, we have that 



e -®(q)x TT (q) 



i \ r°° r°° 

Ug> g) (dx) = J J e -*W <e P*W(fT(t) € dx,L-\t) £ ds)dt 

f'OC f'OQ 

= / / e -^ x e- qs+ ^ x P(H(t) € dx,L-\t) € ds)dt = C/ (g) (dx). 
Jo Jo 

We follow now Bertoin and Doney [5] (see also Kyprianou [20, Th. 7.6 on p. 185]). From Lemma 
6.3 we have 

e^P(X(e q ) >x)= K(q,0) f°° e'^-^U^Jdy) = «(g,0) [°° e'^U^Jx + dy). 

J x JO 

From Kyprianou [20, Th. 5.4 on p. 114] it follows that uj$*(dy) has a nonlattice support. From the 
key renewal theorem (see Kyprianou [20, Cor. 5.3 on p. 114]) the measure U$( q )(x + dy) converges 
weakly to the Lebesgue measure E <s>( q ] H ^ dy (see Kyprianou [20, Th. 7.6 on p. 185]). Thus 

lim e*M*P(X(e q ) > x) - ^ ° ] 



^(q)E iS '(i)H(l)' 



14 



Observe that 



/*oo poo roc 

E^H{\)= te- t E^H(l)dt= e- t E^H{t)dt= xU^(dx) 
Jo Jo Jo 



roo roc roc roc 

/ e-'dt / xP^ q) (H(t) edx)= / e~* / xe^ q >- qs P(H(t) € dx,L~ l {t) e ds)dt 
Jo Jo Jo Jo 



roc roo 

/ e - t EH(t)e*M H ®-* L ~ 1 ®dt= / te 
Jo Jo 



8j3 



From the Wiener-Hopf factorisation (see Kyprianou [20, p. 167]) it follows that 

q-il>(0)=K(q,-0)K(q,0). 

From the convexity of the Laplace exponents <j> and ip we have that k(q,<&(q)) = 2<£(g) > and 
hence n(q, —&(q)) = 0. Finally, 



H(l) = dK ^ 0) 
v ' dp 



/9=-*(?) 



Note that by (6.7) and (6.9), P(X(e q ) > x) = o(e~^ x ) and P(X(e q ) + F(U) > x) = o( e -^ x ). 
Inequalities (6.2) and (6.3) in Lemma 6.1 complete the proof. □ 
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